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The concept of a mean first passage time is used to study the time lapse over which a fissioning 
system may emit light particles. The infiuence of the "transient" and "saddle to scission times" 
on this emission are critically examined. It is argued that within the limits of Kramers' picture of 
fission no enhancement over that given by his rate formula need to be considered. 
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Introduction — Fission at finite thermal excitation is 
characterized by the evaporation of light particles and 
7's. Any description of such a process must rely on sta- 
tistical concepts, both with respect to fission itself as 
well as with respect to particle emission. For decades 
it has been customary to describe experiments in terms 
of particle 1] and fission widths, where the former, Fn, 
is identified through the evaporation rate and the latter 
Ff is given by the Bohr- Wheeler formula Ff = Fbw for 
the fission rate. Often in the literature this is referred 
to as the "statistical model". It was only in the 80 's 
that discrepancies of this procedure with experimental 
evidence was encountered: Sizably more neutrons were 
seen to accompany fission events than given by the ratio 
Tn/rBW (for a review see e.g. 0). A possible enhance- 
ment of that ratio is found if the fission width Fbw is 
replaced by the Fk of Kramers 0. In this seminal pa- 
per he pointed to the deficiency of the picture of Bohr 
and Wheeler in that it discards the influence of couplings 
of the fission mode to the nucleonic degrees of freedom. 
Such couplings will in general reduce the flux across the 
barrier, mainly because of the reduction of the energy in 
the fission degree of freedom Q which may then fall below 
the barrier. This Q is meant to represent the most likely 
path in a multidimensional landscape of shape degrees of 
freedom. 

In Kramers' picture this effect is realized through the 
presence of frictional and fluctuating forces (intimately 
connected to each other by the fluctuation dissipation 
theorem). Presently it is understood that Kramers' "high 
viscosity limit" applies (for a microscopic justification see 
0,13), in which case the rate formula writes 



1 + 7?b - »7b 
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Here, T and E\, stand for temperature and barrier height, 
zua for the frequency of the motion around the minimum 
at Q = Qa and ?7b = (7/(2Mtn))b for the dissipation 
strength at the barrier (at Q ~ Qb) with 7 being the 
friction coefficient and M the inertia. For the sake of sim- 
plicity we will assume these coefficients not to vary along 
the fission path; otherwise the formula must be modi- 
fied 0. For vanishing dissipation strength (Q) reduces 
to the Bohr- Wheeler formula (simplified to the case that 



the equilibrium of the nucleons can be parameterized by 
a temperature). 

Commonly, formula ^ is derived (see e.g.!^) in a 
time dependent picture solving the underlying Fokker- 
Planck equation for special initial conditions with re- 
spect to the time dependence of the distribution func- 
tion 0. Their choice is intimately related to the picture 
of a compound reaction, in that the decay process is as- 
sumed to be independent of how the compound nucleus is 
produced. The latter in a sense represents a nucleus in a 
quasi-equilibrium such that the previous, pre-equilibrium 
stages need not be considered explicitly. This assumption 
is valid as long as the decay of that system takes longer 
than the equilibration time. To some large extent such a 
situation is indeed given at not too high excitations, as 
then the nucleons may stay inside this nuclear complex 
for a sufficiently long time. However, the circumstances 
are less clear with respect to the collective modes, in 
particular to the fission degree of freedom itself — which 
for large damping probably is among the slowest ones 
present. Whereas the corresponding kinetic momentum 
P = MQ may safely be assumed to equilibrate suffi- 
ciently fast, this may not be so for the coordinate Q. 
Thus, assuming the system to be located initially around 
the supposedly pronounced "ground state" minimum of 
the static energy at Q = Qa, the initial width in Q may 
still be at one's disposal. In its true spirit the compound 
picture would suggest taking the equilibrium value, de- 
termined by the temperature and, in harmonic approxi- 
mation, by the stiffness of the potential. Often, however, 
one starts with a sharp distribution of zero width. In 
any case, the current across the barrier needs some finite 
time to build up. This apparent delay of fission was in- 
terpreted 01 as if there was the additional possibility 
of emitting light particles beyond the measure given by 

Tn/rK > Tn/rBW- 

If besides collective motion also particle emission is 
studied explicitly in a time dependent picture, as done in 
the Langevin approach [l^, [ij , such an effect is included 
automatically. Problems arise, however, if one tries to 
imitate this delay in statistical codes which are in use for 
analyzing experimental results. Such codes apply static 
probabilities derived in time independent reaction theory. 
It is not obvious how this method may be reconciled with 
the picture of fission delay, the "transient effect". In the 
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present note we like to shed some light on this problem 
by exploiting the concept of a mean first passage time 
(MFPT). Before we shall come to that we want to ex- 
amine a little closer the time dependent case. We will 
concentrate on over-damped motion, as in this case the 
MFPT can be evaluated from an analytic formula. More- 
over, for slow motion the transient time gets larger, such 
that the feature we want to discuss becomes even more 
obvious. 

Time dependent current across the barrier — In the 
time dependent picture just described the boundary con- 
ditions in Q (and P if present) are chosen to make sure 
that the distribution vanishes at infinity. Calculations of 
the current j{t) across the barrier then typically imply a 
behavior as exhibited in Fig^ In all cases the asymptotic 
value of jh{t) is seen to follow the law Fk exp(— Fk 
shown by the fully drawn straight line. The differences 
at short times are due to the following different initial 
conditions: 
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FIG. 1: The current across the barrier for different initial 
conditions, see text. 

(i) For the dashed and dotted curves the system starts 
out of equilibrium in Q; the dashed curve corresponds 
to the current at the barrier jh{t) = j{Qh,t) and the 
dotted one to that in the scission region jsc(0 = JiQsc, t), 
beyond which the fragments separate. The equilibrium 
is defined by the oscillator potential by which the V{Q) 
around Qa may be approximated. 

(ii) For the fully drawn line the system starts at Qa 
sharp. The obvious delay by about 5 — 10 • 10~^^ sec 
is essentially due to the relaxation of Q to the quasi- 
equilibrium in the well. This feature is demonstrated on 
the right by the i-dependence of the width in Q (exhibited 
in terms of fluctuations of the potential energy CY,qq = 

CiiQ') - {Q?)/2). 

The figure clearly demonstrates remarkable uncertain- 
ties in the very concept of the " transient effect" . First 
of all it is seen, that the "transient" time Ttians, defined 
as the time the current jh{t) needs to reach its asymp- 
totic behavior, depends strongly on the initial conditions. 
Moreover, there is considerable arbitrariness in choosing 
time zero: If the calculation were repeated at some later 
time > Ttrans, the same features would be seen! In the 
end this is due to the very fact that the whole effect only 
comes about because in the initial distribution there are 
favorable parts for which it is easiest to reach the barrier. 



This is demonstrated in Fig|2] There, those points of the 
initial equilibrium are sampled which cross the barrier 
after some given time . On the right a sufficiently large 
Tg was chosen such that greater parts of the initial distri- 
bution have "fissioned". As exhibited on the left, for the 
much shorter time ~ Ttrans, only a small fraction of 
points have succeeded in doing this, namely those which 
started close to the barrier (for under-damped motion 
also more favorite initial momenta would play a role, see 
[T^ IT^l. The vast majority of particles is still waiting 
to complete the same motion but at later times! This 
aspect is important, not only for an understanding of the 
essentials of the concept of the MFPT [Tl[Tl[l6l|. but 
also in respect to the evaporation of neutrons. Indeed, 
even for tk ^ i S> Ttrans there is ample time for them to 
be emitted from inside the barrier. 
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FIG. 2: Samples of initial points which overcome the saddle 
within a time r^: Left part: ~ rtrans; right part Ts — tk- 

The calculations have been performed by simulating 
the Langevin equations exploiting a locally harmonic ap- 
proximation similar to that of 0, |^ for Kramers' equa- 
tion, for the following parameters: T = 3 MeV, Ef, = 8 
MeV, hwa = 1 MeV and r]a = 5. The potential was con- 
structed from two oscillators, one upright and one upside 
down, joined with a smooth first derivative. 

The mean first passage time — Within a Langevin ap- 
proach the concept of MFPT may be described as follows. 
Suppose that at t = particles start at the potential min- 
imum Qa. Because of the fiuctuating force there will be 
trajectories i which pass a certain exit point Qox first 
at some time ti, the first passage time. The mean-FPT 
Tmfpt(Qa Qcx) is defined by the average {ti) over all 
possibilities. In order to really obtain the mean first pas- 
sage time the i has to be removed from the ensemble once 
it has exited the interval at Qcx^ the "particle" can be 
said to be absorbed at Qcx (such that one may speak of 
an "absorbing barrier"). As the potential V{Q) is as- 
sumed to rise to infinity for Q —^ —oo, any motion to the 
far left will bounce back: the region Q ^ —oo acts as a 
" refiecting barrier" . A calculation of the MFPT with the 
Langevin equation is shown in Fig|31by the dashed dou- 
ble dotted curve and seen to be very close to the result 
obtained by exploiting special solutions of the Smolu- 
chowski equation, which we want to address now. 

Fortunately, the Smoluchowski approach allows one to 
derive an analytic formula for the Tmfpt [11,111111 As 
one knows, the Smoluchowski equation represents that 
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of Kramers for over-damped motion. For its solution 
K{Q,t\Qg^,0) the initial condition for the particles to 
start at Qa is given by limt^o K{Q,t\Q a, 0) = 6{Q-Q^), 
which is identical to the one used for the fully drawn line 
of Fig^ For constant friction and temperature one gets 



Tmfpt 



du exp 



Vju) 
T 



dv exp 



V{v) 



T 

(2) 

This expression may be derived as follows: The prob- 
ability of finding at time t the particle still inside the 
interval {~oo, Qox) is given by M^(Qa,i) = /_^^ dQx 
K{Q, 1 1 Qa, 0). Hence, the probability for it to leave the 
region during the time lapse from t to t+dt is determined 
by —dW = —{dW{Qa.,t)/dt)dt, such that the average 
time becomes Tnifpt(Qa Qcx) = ~ J tdW which turns 
into 



'^mfpt {Qa ^ Qc 



K{Q,t\Q,,0)dQdt 



dUj(Qex,t|Oa,0). 



(3) 



These formulas are associated to the special boundary 
conditions with respect to the coordinate mentioned be- 
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FIG. 3: MFPT for a cubic potential normalized to its asymp- 
totic value, shown by the solid, dashed, dotted-dashed and 
dotted curves which correspond to T/Eb — 0.1,0.2,0.5,1.0. 
The dashed double dotted curve represents a calculation 
within the Langevin approach for T / Ef, = 1. 

fore, the reflecting barrier at Q — > — oo and an absorbing 
barrier at Qox- In particular for the latter feature it is 
not permitted to use in © the currents shown in Fig^ 
Inserting them blindly would indeed lead to expressions 
for Tnitpt in which the Ttrans appears jT^. This is in 
dear distinction to the correct form Actually, the 

derivation of ^ involves proper solutions of that equa- 
tion which is "adjoint" to the Smoluchowski equation, 
and which describes motion backward in time. In FigO 
we show the dependence of Tinfpt(Qa — > Qox) on Qox as 
given by Q calculated for a cubic potential. Evidently, 
the MFPT needed to reach the saddle at Qb is exactly 
half the asymptotic value. The latter may be identified as 



the mean fission life time T{ = Tinfpt(Qa — > Qox ^ Qb)- 
For the typical conditions under which Kramers' rate for- 
mula is valid for overdamped motion, the identity of 
Tf = tk to the asymptotic value of the MFPT can be 
proven analytically |l5| . Another remarkable feature seen 
in FigOis the insensitivity of the MFPT to the exit point 
for small and large Qcx- Actually, in clear distinction to 
the transient time the MFPT is also insensitive to the 
starting point. This latter property shall be exhibited in 
a forthcoming paper |l3j |. 
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FIG. 4: The MFPT for the case of small and vanishing bar- 
riers for T = 2, 3, 4 MeV from top to bottom. The Qsc corre- 
sponds to a scission point 20 MeV below the barrier. 

As a most interesting feature, the MFPT can be cal- 
culated also for cases of small barriers where Kramers' 
formula does not apply. We show in Fig0]results of eval- 
uations of formula ||2J) for E\j = 1 MeV as well as for a 
practically vanishing barrier. It is seen that even in the 
latter case the Tinfpt(Qa Qcx) reaches a plateau for suf- 
ficiently large Qox- This asymptotic value, however, is no 
longer determined by tk- Nevertheless, the Tmfpt seems 
to be long enough for neutrons to be evaporated before 
scission. This may be seen as follows. The neutron width 
typically is of the order of 1 ... 3 MeV (for medium heavy 
nuclei 10]). According to the 7/C roughly increases 
linearly in T, being about 1 h/MeV at T = 2 MeV and 
about 3 h/MeV a.tT — 3 MeV. Taking the values of Tmfpt 
from the right part of Fig^and close to the plateau one 
gets a width Fmfpt — ft/Tmfpt of the order of 0.9 ... 2.3 
MeV, and thus comparable to the neutron width. Fig^ 
also shows that a small increase of the barrier by 1 MeV 
enlarges the Tmfpt drastically. Discussion — It should 
be evident from the previous discussion that in the very 
concept of the MFPT there is no room for a transient 
effect. After all, formula ^ is based on exact solutions 
of the transport equation which satisfy the same initial 
condition as those used for the plots in Fig^ — albeit 
different boundary conditions in coordinate space. More- 
over, as exhibited in Fig[5]the evaluation of the MFPT 
takes into account an average over all initial points, as 
is warranted by the definition of the MFPT through the 
probability distribution —dW. Contrasting this feature, 
and as outlined in the second section, the transient effect 
only represents a minor part of the initial population, 
namely that one which reaches the barrier first. Discard- 
ing the rest implies ignoring the many particles which are 



4 



still moving inside the barrier for times typically much 
longer than Tti-ans- Hence, neutrons from deformations 
corresponding to that region may not only be emitted 
within Ttians but within TmfptiQa ~* Qh), which turns 
out to be just half of the total fission time tk- Of course, 
this discussion shows that it is also not correct to argue 
in favor of "additional" neutrons which might be emit- 
ted within the saddle to scission time Tgsc introduced in 
[l8| . As one may guess from Fig|31 hke the Ttrans, the 
Tgsc does not appear to be in accord with the MFPT ei- 
ther: The time the fissioning system stays together is not 
determined by motion in the immediate neighborhood of 
the barrier. On average it takes half the full decay time 
to move beyond the Qb to the Qex at which the Tmfpt 
reaches its plateau value. 

These findings suggest that one simply estimates the 
emission rate of neutrons over fission from the ratio 
Tn/rK — provided one may trust the potential to be 
of the simple form underlying the rate formula . Any- 
thing else does not seem to be in accord with an appro- 
priate application of Kramers' or Smoluchowski's equa- 
tions. This does not rule out other, complementary ef- 
fects which originate in more complicated situations. For 
instance, in case that in the scission region the poten- 
tial becomes flat again or even develops a minimum the 
system is forced to stay there longer than given by the 
Tk of eq.(^ — implying additional time for evaporat- 
ing neutrons. Likewise it is conceivable that the initial 
stage of the whole reaction is to be described with a dif- 



ferent transport model. Such modifications are already 
suggested when the average neutron emission time 
becomes comparable to or even smaller than the relax- 
ation time Tijiicro for the nucleonic degrees of freedom as 
a whole. Transport equations are justified only if this 
Tmicro is the Smallest time scale present, in comparison 
to both neutron emission as well as to collective motion. 
The Tniicro tums out to be of the order of f — 2 • 10^^^ sec, 
no matter whether it is estimated within linear response 
theory with coUisional damping or within a random ma- 
trix approach (see 0). Applying the Weisskopf estimate 
for the neutron emission time (or modified versions of 
it) [3 ^-t large temperatures one easily gets values of 
Tn of the order of or smaller than Tmicro- This reflects 
a situation of pre-equilibrium rather than that assumed 
in the quasi-static picture necessary for the application 
of Fokker-Planck equations. In conclusion we may say 
that deviations of experimental results from the stan- 
dard value of Fn/FK ought perhaps to be understood as 
a strong indication of the relevance of these complemen- 
tary effects, which unfortunately have for the most part 
been unconsidered. This might require one to re-examine 
analyzes of experiments which over the past decade or so 
have followed the conventional line. 
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